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I. INTRODUCTION 

The loss of quantum coherence due to unavoidable in- 
teractions of quantum systems with the surrounding en- 
vironment is known as decoherence. It represents the 
main obstacle to quantum computing and quantum infor- 
mation processing^'^'^. The environment destroys quan- 
tum interferences of the central system within time scales 
much shorter than those typically characterizing dissipa- 
tion''. The unwanted effect of decoherence reduces the 
advantages of quantum computing methods by produc- 
ing errors in their outcomes. Different strategies, such as 
error-correcting codes, are adopted to overcome this dif- 
ficulty^i^ii^. Great scientific effort has been devoted to 
the understanding of the process of decoherence in quan- 
tum systems, mainly focused on solid state spin nanos- 
tructures. These systems seem to be the most promising 
candidates that can be efficiently used in quantum infor- 
mation processing and comput at ion^ d ° d ^ . 

Several models were proposed to study decoherence 
of single and multi-spin systems interacting with a sur- 
rounding environmen1>i^. Very often, the derivation of 
the reduced dynamics involves complications and diffi- 
culties that can be overcome in many cases by making 
recourse to approximation techniques. In particular, the 
Markovian approximation together with the master equa- 
tion approach turns out to be very usefu l '^'^d^ . However, 
any approximation method is inevitably based on some 
assumptions which do not necessarily reflect the actual 
properties of the composite system. Moreover, many re- 
alistic spin systems exhibit non-Markovian behavior for 
which the standard derivation of the master equation 
ceases to be applicable. The non-Markovian dynamics 
of a central spin-system coupled to a spin environment 
has been investigated by many author a^^d^d'^d^d^ . 

In general, the course of the decoherence process de- 
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pends on the intrinsic properties of the bath such as 
temperature, polarizations, and quantum fluctuations. 
At low environmental temperatures, the dominant effect 
arises from the contributions of localized modes such as 
nuclear spins^". In quantum dots, the decoherence of 
the central spins is mainly caused by the hyperfine cou- 
pling with the surrounding nuclear spins. The effect of 
bath polarizations and external magnetic fields on the 
decoherence of electron spins in quantum dots has been 
investigated by Zhang et al—. 

In the following paper we study the dynamics of a spin- 
^ particle interacting with a large spin environment in 
thermal equilibrium. In Sec. II we introduce the model 
Hamiltonian together with the initial states of the central 
spin and the environment. In Sec. Ill, we calculate the 
exact time evolution operator of the composite system 
and we derive the reduced density matrix of the central 
spin. Sec. IV is devoted to the case of an infinite number 
of spins in the bath. We study the long-time behavior 
as well as the short-time behavior of the reduced density 
matrix, and we discuss the effect of the magnetic field and 
the bath temperature on decoherence. A short conclusion 
ends the paper. 



II. THE MODEL 

We consider a central spin-^ particle coupled to a spin 
bath composed of A'^ interacting spin-i particles in ther- 
mal equilibrium at temperature T. The spin operators 
corresponding to the central system are denoted by S'^ 
with i — x,y,z, those associated with the bath con- 
stituents are denoted by where k = 1,2,...,A^ and 
i — x,y, z. We assume that the central system as well as 
every spin in the bath couples to all other spins through 
long-range anisotropic Heiscnberg interactions. More- 
over, an external magnetic field of controlled strength fi is 
locally applied to the central spin along the z direction. 
Under the above assumptions, the model Hamiltonian 
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can be written as 



Sl.±iSy as follows 



H = Hs 



H. 
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where Hs and Hb are, respectively, the Hamiltonian op- 
erators of the central spin and the surrounding environ- 
ment. The coupling between the open system and the 
bath is described by the Hamiltonian Hgg. Explicitly, 
we have 



Hs = 2^5°, 



(2) 
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where 7 and a are the coupling constants of the central 
spin to the environment, g stands for the strength of in- 
teractions of spins in the bath, and A is the anisotropy 
constant. The coefficient 2 in front of /x, 7 and a in 
Eqs. ([2]) and ([3]) is introduced for later convenience. 
Furthermore, we have rescaled the above interaction 
strengths with appropriate powers of the number of spins 
in the environment in order to ensure good thermo- 
dynamical behavior, namely an extensive free energy. 
Obviously, a more realistic model would include site- 
dependent interactions. 

Note that in the case where 7 = 0, Hsu reduces to 
Heisenberg XY Hamiltonian which was recently used to 
model the coupling of one and two qubits to star-like 
environment^i^ii^i^i^. Moreover, when 7 = 0; we sim- 

_^ N ^ 

ply have Hs b — -f^ ^ 5"% which should be compared 

^ i=l 

with the Hamiltonian of the hyperfine contact coupling 
of electron spin to the nuclear spins in quantum dot . In 
Ref. m, the Hamiltonian hs = T.i>j QijiS^S^ - iSlSi) 
was used to model the intra-bath dipolar coupling be- 
tween nuclear spins in quantum dot. If we assume uni- 
form coupling between nuclear spins, i.e. all the gij 
are the same, then the operator Kb (with rescaled cou- 
pling constant) becomes equivalent to Hb in the case 
where A = — 2. It should also be noted that the bath 
Hamiltonian Hb is very close to that of the isotropic 
Lipkin-Meshkov-Glick model2ii2^. There, the magnetic 
field globally applied to all spins plays the role of the 
anisotropy present in our model. This can be better seen 
by applying mean field approximation to the longitudinal 
term of Hb- 

The Hamiltonian operators Hb and Hsb can be 
rewritten in terms of the lowering and raising operators 
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HsB^^siy^si + ^\siy^si + s'_y^s\ 
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By introducing the total angular momentum of the bath 

N _^ 

J — ^ S"*, together with the corresponding lowering 

1=1 

and raising operators J±, it is possible to put the above 
Hamiltonian operators into the following form 



Hsb = 
Hb = 



27 



N 
9 

2N 



S^Jz 



K + 2A J: 



a 
N 

2 (2 + A)7V 



(7) 
(8) 



Here, is the z-component of the total angular mo- 
mentum J, and we have introduced the operator K — 
J+J^ + J^J+. From here on, we shall neglect the con- 
stant (2 + A)g/4: appearing in the expression of Hb since 
it has no effect on the dynamics of the system. This can 
be done by redefining the energy origin of the spectrum 
of the bath Hamiltonian. 

The spin spaces corresponding to the central spin and 
the environment are given by and (C^)*^^, respec- 
tively. The latter space can be decomposed as a direct 
sum of subspaces C''^ each of which has a dimension 
equal to dj = 2j + 1 where < j < we take N 



even), namely (C 



u{N,j)C'^'. The degener- 



acy v{N,i) is given by= 
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It is worth noting that the bath Hamiltonian can be ex- 
pressed in terms of the operators and Jj as Hb = 
;^[J2 + (A - 1)J2]. Therefore, the operator Hb is di- 
agonal in the standard basis of (C^)*^^ formed by the 
common eigenvectors of and Jz which we denote by 
|j, m) where —j < m < j. In this basis, the eigenvalues 
of the operator K are simply given by 2{j(j -I- 1) — m^) 
(we set h = 1). 



III. REDUCED DYNAMICS OF THE CENTRAL 

SPIN 



In this section we derive the exact time evolution of the 
central spin for finite number of environmental spins. As 
usual, we introduce the time evolution operator U(i) = 
together with the total density matrix operator of 



-iHt 



the spin-bath system, ptot{t). The initial value of the 
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latter is denoted by ptot (0) . The evolution in time of the 
composite system is unitary, its density matrix at any 
moment of time is given by 



PtotW = u(Optot(o)ut(t). 



(10) 



The reduced density matrix of the central spin can be 
calculated by tracing ptot (t) with respect to the environ- 
mental degrees of freedom, namely 

pit) = tTBlptotit)}. (11) 

This can be explicitly written in terms of bath states as 
pW -I]K^,j)(j,"^|Ptot(i)li,"i)- (12) 

In order to solve the time evolution problem (fTO|) . one 
needs to calculate the exact analytical form of U(t) and 
to specify the initial density matrix. 

A. Initial conditions 

Initially, the central spin is assumed to be uncorrelated 
with the environment. The corresponding total density 
matrix is given by the direct product /Otot(O) — p{0) (g) ps 
where p{0) and pB are, respectively, the initial density 
matrices of the central spin and the bath. In the standard 
basis composed of the eigenvectors |— ) and 1+) of the 
operator S^, Ps(0) takes the general form 



(13) 



where pn and p22 are positive real numbers which satisfy 
pii + P22 — 1. For instance, if at < = the central system 
was in the state 



\m) = a\-)+b\+), 



(14) 



where a and b are complex numbers satisfying |ap + |6p = 
1, then p^i = l«P and p'(2 — ab* . 

Alternatively, p{0) can be expressed in terms of the com- 
ponents of the Bloch vector A = (Ai, A2, A3) as 

„(n)-U 1-^3(0) Ai(0)-^A2(0)\ 
2 U(0) + *^2(0) 1 + A3(0) )' 

with the condition |A| < 1; the equality holds for pure 
initial states. We shall use both representations of the 
density matrix throughout the paper. 

At t = 0, the spin bath is taken in thermal equilibrium 
at finite temperature T. Its density matrix is given by 
the Boltzmann distribution 



PB 



-PHb 

Zn 



(16) 



where /3 = 1/T (we set ks = 1), and Zn = tr^e-^^^ is 
the partition function of the bath. Clearly, ps is diagonal 
in the standard basis m)} from which it follows that^i 

(j,m|pB|j,m) = (17) 



and 



Zn = Y,v{N,j) e-#b-0+i)+(A-i)m=]_ (^g) 



In the case of the isotropic Heisenberg model, i.e. when 
A = 1, the above expression simplifies to 



ZN=Y.y{N,3){2j + l) e 



(19) 



In the extreme case of an infinite temperature (/3 — >■ 0), 
the density matrix of the bath reads 



Pb{T = 00)^^, 



(20) 



which corresponds to a completely unpolarized spin bath. 
In the previous expression 1 b stands for the unity matrix 
in the bath space. 



B. Time evolution operator 

Let Uij denote the components of the time evolution 
operator U in the basis {|— ), corresponding the cen- 
tral system space. Therefore, we can write 



U|-) = C/ii|-) + C/2i|+), 

U|+> =C/l2|-)+t/22|+). 



(21) 
(22) 



On the other hand, the operator U satisfies the 
Schrodinger equation 



z-U|±)=HU|±) 



(23) 



Substituting Eq. (|2Tj) into Eq. (|23l) yields the following 
system of coupled differential equations 



■TT ^-^^ TT 

1U21 = -^=Uii 



N 



Similarly, from Eq. (gH) and Eq. ^ we obtain 



f/21, 
N 

(24) 

C/21. 
(25) 
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U2: 
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(27) 
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Since U(0) = I2 (Xi 1_b, one gets the initial conditions 

C/ii(0) = C/22(0) = li3, C/i2(0) = C/2i(0)=0. (28) 

The difficulty with solving the above set of differential 
equations resides in the fact that the coefficients of the 
operator variables t72i and U12 are not diagonal and do 
not commute with those of Uu and 1/22- Nevertheless, 
as we shall see bellow, this problem can be overcome 
by transforming these equations into new ones involving 
commuting diagonal operators. Indeed, by making use of 
the change of variables (see Ref . il^ for a similar method) 



Vn' ' 2JV 
1-1: 



iK+2Aj't)]tf 



C/21 = J_c-'[-(^+^) + ^(^+^^^^')l*C/2i, 

U22 = e-'t(''+w)+^(^+2'^^--)l*C/22, 
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^12, 



(29) 
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and taking into account the commutation relations 



[J„ J±] = ±J±, [Jl J±] = ±J±(2J, ± 1) 



and 



[if,J±] = T2J±(2J,±1), 



(33) 



we obtain 



iUn = 
iU 21 = 
iU22 = 
1U12 = 
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J+J-U21, 
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J-J+U12, 
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(1-A))(J.4) 



(1-A) J, 



U21, 



U12 
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Now, the coefficients in front of the new operator vari- 
ables Uij are diagonal in the common eigenbasis of 
and Jz, whence the standard method of solving systems 
of differential equations can be easily applied. Combin- 
ing the above relations leads to the following second order 
homogeneous differential equations for the operators U21 
and U12 



U21 + 2i 



vn n"- 'jy 2 



U21 



-f ^J+J_C/21=0, 
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+ —J-J+Uu 
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N ■ N 
= 0, 



(1-A) J, 



Ui 



which admit the following solutions 



C/21 = 2iCie- 



"^*sin( 



U12 = 2iC2e"'^^* sin( ty^M2 



i(i\/M2). 



(38) 
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(39) 

(40) 
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Here, Ci and C2 are some diagonal operators to be de- 
termined and we have 
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N N 
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(1-A))(J.4), 



M, = F,^ + -J+J. 



F2 



M2 = F4 
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-J_J+ 



(42) 

(43) 
(44) 

(45) 



Integrating the right-hand side of Eq. (|40|) gives 

c7n = -2C,e--^*^ 
a 

iFi 



-C3. 

(46) 



The constant operators Ci and C3 can be determined 
using the initial conditions (j28p and the unitarity condi- 
tion for the time evolution operator, which yield Ci = 



5 



tj {2^/NAIi)1b and C3 = 0. Hence, we obtain 



Unit) = e 



-iGit 



(47) 

UMt) = -^J-^^e'^^^^^in{t^), (48) 
where 



r - 7,5 



i^ + 2Aj2)+2(l-A)(j,-i) 

(49) 



Following the same method, we find that 



U22{t) = e-'^^*[cos(tv/M2)+^^sin(tyM2) 



^-W = -J.^e---sin(.yMg), 
where 



(50) 
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Go = 
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(i^ + 2Aj2)+2(A-l)(j, + i 



(52) 



It is easy to see that the operators Gi and G2 are diagonal 
in the common basis of and J,. 



C. Reduced density matrix 

Having determined the exact analytical form of the 
time evolution operator, we are able to calculate the re- 
duced density matrix of the central spin. Indeed, from 
Eqs. (fTO|) and pT|) . and by making use of the trace prop- 
erties of the lowering and raising operators J± , we find 
that 



p02trB(e-^^-C/2*i[/i2 



Pi2(i) 



Zn 



(53) 
(54) 



Furthermore, with the help of the commutation rela- 
tions we can easily prove that J^Fi — ~F2J-, 
and J- Ml — M2J-. Using the latter equalities, one can 
check that the time-dependent components of the Bloch 
vector are given by 



X,it) = -^tr5{^^e-|#[^+2A,/.?+(i-A)(2,7.-i)] ,i^2 (^t^y^^^l^^^ _ ^)(2J, - 1)] } 



2 



+ A3 



NMi 

(0)|l - AtrB|4^e-w[^+2Ajf+(i-A)(2j.-i)] sin2(t^)cosh 



|^(1-A)(2J,-1) 



Ai(t) = trB|(^Ai(0)cos(rJi) -I- A2(0)sin(m))A- (^Ai(O) sin(f7i) - A2(0) cos(m))s|, 
|(Ai(0)sin(m) - A2(0)cos(rJt))A-f (Ai(0) cos(f)t) -t- A2(0) sin(m))B|, 



where 



-tr_B 



f7 = ^(A-l)J„ 



A = ■! e 2JV 

Zn 



1 

Zn 



^''1^ ^^ = sin(^iyM^^cos(^tyM2)-;^^=sin(^iyA/2)co |. 
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From here on, the parameters fi and 7 will be given 
in units of the coupling constant a. The behavior of the 
component A2(t) does not significantly differ from the 
one corresponding to Xi{t). Throughout the remainder 
of the paper we shall deal with the latter component and 
restrict ourselves to positive values of the anisotropy con- 
stant A. 



Depending on the nature of interactions within the 
bath, we can distinguish two different cases. The first one 
corresponds to positive values of 5, i.e. antiferromagnetic 
couplings between the constituents of the environment. 
In this the number of spins increases, the plots 

saturate and a nontrivial limit exists as shown in Fig. [T] 
This will be investigated in the following section. The 
other case corresponds to negative values of g, i.e. fer- 
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FIG. 1: Time evolution of the components X3{t) and Xi{t) for different values of the number of spins in the environment: 
N — 100 (dotted lines), TV = 200 (dashed lines), and A'^ — 400 (solid lines). The other parameters are 7 = 0, (? = l,/3 = 0.5, 
A = 0, and fj, = a. The initial conditions are A3(0) — ^, Ai,2(0) — |. 




«t at at at 

(a) (b) 

FIG. 2: The Gaussian decay of the Bloch vector components X^it) and Xi{t) in the case of ferromagnetic interactions: (a) 
N=100 and (b) N=200. The plot on the left of each subfigure corresponds to A3(t), the one on the right corresponds to Xi{t). 
The other parameters are 7 = 2a, [3 — 1, g = —5, A — 0.5, ~ 0, A3(0) — ^, and Ai,2(0) — |. 



romagnetic couplings within the bath. When A < 1 the 
components of the Bloch vector exhibit in general Gaus- 
sian decay accompanied by fast damped oscillations even 
when the strength of the magnetic field is very weak. In 
contrast to Xi{t), the component X^it) decays faster as 
the number of spins increases. When the latter is small, 
Xslt) may revive to decay again and so forth. The nu- 
merical simulation shows that the details of the time evo- 
lution of the reduced density matrix are rather complex 
and depend on the different values of the parameters of 
the model, including the number of bath spins. For ex- 
ample if we set A = 0, we observe that the oscillations are 
quickly suppressed with the increase of the strength of the 
magnetic field, or the value of the coupling constant 7. 
In this case, the components A2,3(t) do not vanish at long 
time scales; the corresponding asymptotic values depend, 
however, on N in contrast to the antiferromagnetic case. 
For large values of the coupling constant g, the compo- 
nent Ai (t) quickly decays whereas A3 (t) oscillates around 
zero with large amplitudes (typically of the same order of 
magnitude as the corresponding initial value). We also 
notice that the frequencies of the damped oscillations in- 
crease with the increase of the number of bath spins as 
shown in Fig. [21 Roughly speaking, when A > 1, the 
behavior of the components of the Bloch vector is quiet 
similar to the antiferromagnetic counterpart. For exam- 
ple when 7 = 0, the components Xi(t) show saturation 
behavior with respect to the number of spins A'^; their 



asymptotic values are different from zero. 

In order to explain the differences between the behav- 
ior of the reduced density matrix in the ferromagnetic 
and the antiferromagnetic environments, we note that in 
the latter case, the form of interactions favors antiparallel 
spins. This is the reason for which the ground state of the 
antiferromagnetic bath, |5'g), is equal to |0, 0). On the 
contrary, ferromagnetic interactions force the spins in the 
bath to align along an arbitrary direction in the space. 
In this case belongs to the subspace C^"*"^ spanned 
by the state vectors | ^ j corresponding to j = ^ ■ For 
instance, when A > 1, the ground state of the bath turns 
out to be doubly degenerate, namely |\E'g) = Iyj^t)- 
For A < 1, we simply have Iv^g) = It'*^)- However, 
when A = 1, the ground energy of the bath is indepen- 
dent of the quantum number m; the degeneracy of |^'g) is 
equal to -I- 1. Hence we conclude that the Hamiltonian 
Hb displays quantum phase transition at A = 1 . This is 
the reason for which the reduced dynamics depends on 
whether the anisotropy constant is less or greater than 
one. Note that the mean value of is close to zero 
in the case of antiferromagnetic interactions within the 
spin bath in contrast with the ferromagnetic case where 
(J^) ^ N^. Obviously, the central spin decoheres less 
if the spin bath, to which it couples, is characterized by 
a total angular momentum close to zero. At zero tem- 
perature, the antiferromagnetic bath occupies its ground 
state 1 0,0) which is an eigenvector of Hb, and satisfies 
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HsbI^) ® |0, 0) = 0. Hence, the central spin remains de- 
coupled from the bath if the initial state factorizes: the 
two-level system preserves its coherence regardless of the 
number of environmental spins. Let us now consider the 
case where 7 = and A > 1. At low temperatures, the 
total angular momentum of the ferromagnetic bath has 
the tendency to be directed along the z direction. Since 
the central spin couples to the bath through Heisenberg 
XY interactions (7 = 0), we end up with a situation 
quiet similar to that where g > 0. The above results 
show that properties of the bath at zero temperature af- 
fect the behavior of the reduced dynamics when T > 0. 
At infinite temperature, the ferromagnetic and antiferro- 
magnetic environments become completely unpolarized; 
the reduced dynamics displays the same behavior in both 
systems as N increases. 



IV. THE LIMIT N 



which is a special case of / ^zn+ig-at _ 20"+^ ' "^here 



n = 0, 1, 2, and the real part of a satisfies Re(a) > 0. 

The operator J^/Vn also converges to a real random 
variable m (to be differentiated from the eigenvalue m) 
when N 00, with the probability density function 



(65) 



m t-^ \l — e 

TT 



For example, consider the operator e 'vw ^ and let us 
calculate 



N 



tr e n/Tv ^ . 



(66) 



k=l 



The trace under the product in the right-hand side of 
the above equation can be easily evaluated as 2 cos(-^). 
Consequently, 



This section is devoted to the case of an infinite num- 
ber of spins in the environment, i.e. the case N — > 00. 
We investigate the effect of the bath temperature, the ex- 
ternal magnetic field, and the anisotropy constant on the 
reduced density matrix of the central spin. To this end it 
should be noted that the trace of the operators J± / \/N 
together with J^/VN is identically zero, namely 



0. 



(61) 



A more general property of the trace of the lowering and 
raising operators can be expressed as 



lim 2-"trB , . , , 

i=l 



'{n 



} 



N 



hm 2--tr4TTf4^)"f4^)"|-?' 



(62) 



tr.{e-^-}=2^[cos(^) 



N 



(67) 



Expanding the cosine function in a Taylor series and tak- 
ing the limit N ^ 00 yield 



lim 2"^tri3 

N^oc 



lim 



1 



72^2 

2N 



N 



On the other hand we have 



00 



2.2 



TT 

00 



(68) 



(69) 



which is in agreement with Eq. 
can infer that 



In particular we 



JV 



r(7i 4- i) 



lim 2-''ti-BiJ./VNy "^ ^Z 'f^' ' (70) 



2"V7r 



where n — n.i is positive integer; the trace vanishes for 
1=1 

all the cases in which J+ and J- appear with different ex- 
ponents. This means that J±/\/N are well-behaved fluc- 
tuation operators with respect to the tracial state. Hence 
in the limit N ^ 00, the operator Jj^/^/N converges to a 
complex random variable z with the probability density 
functioni^ 



(63) 



where r(z) is Euler gamma function. We shall use the 
latter results when we investigate the short-time behav- 
ior of the reduced density matrix in the case where 7 is 
different from zero. 

One can check that for large values of N , 



tr, 



Here, we wish to mention the similarity that exists be- 
tween relation (l62t and 



00 

4 y" t dtt 





2" 



(64) 



For odd powers of Jz, the left-hand side of the above 
relation vanishes as N increases; the right-hand side is 
always zero. Eq. (|7ip simply implies that the opera- 
tors J±J^/N and J^/VN become uncorrelated under 
the tracial state at large values of N . Note that the 
above state corresponds to a bath of N independent 
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spin-i particles, i.e. the state of maximum entropy. 
In the hmit of large number of spins such a bath has 
the tendency to behave as a classical stochastic system. 
The scaled bath operators Jq/a/ZV (where a = x,y,z) 
converge to independent commuting random variables. 
For instance, we can easily show that the trace over 
the environmental degrees of freedom of the operator 



exp 



ct 

Vlv 



aiJx + a2Jy + a^J; 
e — \/±T, is given by 2^ | cosh 



where ai 



2,3 



and 

N 



If we expand the cosh function in Taylor series and take 
the limit N — > oo, as we did in Eq. ([S5)) . we end up 



with the result exp[^-|-(a^ 



■a\)]. The latter can be 



obtained by multiple integration over three independent 
random variables each of which has the same probability 
density function as m [see Eq. (IMl) ]. It follows that the 
random variables z and m can be treated as independent 
in the limit N ^ oo. 

From the above discussion, we can conclude that 



lim 2"^trB 



J±Jzr Jz 



N 



hi. 



2\3/2 



dm / dzdz*f{\z\'^,m) e^2(™ +NI 



(72) 



at least for bounded functions / : C x M ^ M. The latter 
relation has been numerically checked for large number of 
functions; the agreement between its two sides is perfect. 
In fact, the class of functions for which the integral in 
the right-hand side of Eq. (TtS]) exists contains all the 
fimctions having the form e"^"'^' ^h{\z\'^,m) where 
h is bounded and a and b are complex numbers satisfying 
Re(a) > —2, Re(6) > —2. If the latter conditions are 
not satisfied then the integral does not converge. This 
is the reason for which we shall restrict ourselves to the 
antifcrromagnetic case where g and A are positive. 

Under the above assumptions, it is possible to evaluate 
the quantity 



lim 



N- 



2\3/2 



dm / dzdz*e-(2+9/3A)m^-(2+<,/5)|.|^ 



(73) 



by making use of the polar coordinates (r, 0) where z 
re*"^. A straightforward calculation yields 



2y/2 



{2 + g(3)V2 + gPA (2 + 5/3) 



(74) 



when A ^ 0. Obviously, if gjS = then Z = 1. The 
agreement between the right-hand side and the left-hand 
side of Eqs. (TtS]) is illustrated in TableUwhere we display 
Zn at different values of N and compare it with Z 
for g = 2, A = 5 and /? = 1; the agreement is clearly very 
good for N = 5000. 

Let us now focus on the general structure of Eqs. ([55)1 - 
(157)1 . Clearly, we need to evaluate terms having the gen- 
eral form 



TABLE I: 2J^Zjv at different values of AT for g = 2, A = 5 
and /3 = 1; ^=0.204124. 



N 


10 


100 


1000 


5000 


2-^Zjv 


0.203026 


0.203997 


0.204111 


0.204122 



Zjv""L^V N '^/]v^J 2-^Zn 
As ^ cxD, the previous quantity tends to 



(75) 



1/2/* /* 

(/) = ^"^4(-) / dm / rdrfir\m) e-2(™'+'-') 



(76) 

This is permissible since the functions of interest appear- 
ing in Eqs ([55])- ([57]) fulfil all the conditions mentioned 
above. Note that the factor 4 in Eq. (|75|) appears af- 
ter performing the integration with respect to the polar 
coordinate (f) ( this actually follows from the symmetry 
with respect to the z direction). It is also quiet interest- 
ing to notice that the behavior of the central spin when 
—2 < g(3 < and —2 < gPA < is similar to that where 
g > and A > as indicated by the conditions on the 
convergence of the integral in Eq. (17^ . 



A. The case 7 = 

Let us assume that the coupling constant 7 is equal 
to zero. First of all, it should be noted that, although 
the operator J^/VN converges to a random variable, we 
can neglect the contribution of Jz/N when A^ becomes 
very large. This means that in the limit N 00, the 
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FIG. 3: Evolution in time of X^it) and \i{t) for iV = 400 
(dotted lines), and N —> oo (solid lines); the dashed lines 
correspond to the asymptotic values. Other parameters are 
7 = 0, = 1, /? = 5, A = 0.5, fi = 0.4q, A3(0) = i and 
Ai,2(0) = |. 



quantities Mi. 2 do not depend on the random variable 
to; the sinh (sin) and the cosh (cos) functions appearing 
in Eq. ([55)1 [Eqs. ([55)1 - ([Ff|) ] should be replaced by zero 
and one, respectively. We only need to integrate with 
respect to the random variable z since the integrals with 
respect to m occurring in the numerator and denominator 
of Eq. ([7S)l cancel each other. One then concludes that 
the anisotropy constant A has no effect on the dynamics 
of the central spin when TV — > 00. This is due to the 
fact that HsB simplifies to Heisenberg XY Hamiltonian. 
Only transverse interactions contribute to the reduced 
dynamics when N is sufficiently large because Jz/VN 
and K/N (or equivalently J±J^/N) become practically 
uncorrelated under the tracial state [see Ea. ([7T|) ]. 

Hence, in the limit of an infinite number of spins within 
the bath we obtain 

A3(i) = A3(0)(l-r;(0), (77) 

where 

sin^ (tJu'^ +rA 
m = (2r- y^^2 ^ ,-.m.-.A„.-l). (78) 

Note that the time variable is now given in units of a. 
We show in the appendix that the above function can be 
written as 



ri{t) = 1 — cos( 



erf 



Ai(.9/9 + 2) -it 



-erf 



Ai(.g/3 + 2)+ it 



} 



,[(2+g/3)p2 



VWT2 

{g(3 + 2)fi^ e(9/3+2)A''r[0,(g/3 + 2)/i2]+/i2(g/3 + 2) e(^"^+2)^'Re|r [o, (5/? + 2)/i2 + 2^ii 
H^M{t;fi,P), 



} 



(79) 



where 



erf(z) 



r(a,z) = / t 



"* dt. 



'"-^e'^dt. 



(80) 



(81) 



the appendix. 

The remaining components of the Bloch vector are 
given by 



Mt) = Ai(0) [at) + ^fj{t)] +A2(0)e(t), (82) 
A2(t) = A2(0) [at) + ^rjit)] -Ai(0)e(t), (83) 



are, respectively, the error and the incomplete gamma 
functions^^. The function Ai is given by Eq. (|A.10[) of where 
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VqBX 



2 + .g/9 



erf 



erf 



(84) 



and 



sinf 2^^//i2 + r2 



^V ^(2 + g/3) e[(^+^'^)^^-^l {erf r^^^^ + " --rA'C?/? + 2) + 



—erf 



(85) 



The asymptotic behavior of the reduced density matrix 
can be easily determined as follows. Let us begin with 
the simplest functions namely C(t) and ^(i). Their limits 
when t — s- oo are equal to zero which immediately follows 
from the Riemann-Lebesgue lemma 

lim C(i) = lim i{t) = 0. (86) 

t — >oo t — *oo 

The same lemma can be applied to the function 77 (i) after 
some simplifications of the integrals of interest as shown 
in the appendix. Only one term survives the above ap- 
proach when t goes to infinity, namely 

lim r/(t) = l-r/°°, (87) 

t — 'OO 

where 

= P?{g^ + 2) e'''(f''+2)r(0, p?{gP + 2)) . (88) 

Hence, the asymptotic behavior of the reduced density 
matrix can be expressed as 

lim \(t) ^Xo- T]°°WX{Q) (89) 

t — >-oo 

with 

^ /Ai(0)\ a o\ 

Ao= 2 I A2(0) I , W = I i I . (90) 

The evolution in time of the components A3 (t) and Ai (t) 
is shown in Fig. [3] for N — 400 spins in the environ- 
ment, along with the corresponding infinite case and the 
asymptotic limits obtained in Eq. ([55)1 . We can see that 
the off-diagonal elements of the reduced density matrix 
show partial decoherence. At low temperature, the rele- 
vant bath states are those with low energies (i.e. j close 
to zero). In this case, the central spin is weakly cou- 
pled to the bath and hence preserves most of its coher- 
ence. At high temperature, the two-level system becomes 
more correlated with the bath which, however, behaves 
as a system of independent uncoupled particles. Thus 



quantum fluctuations within the antiferromagnetic spin- 
environment reduce the effect of the decoherence of the 
central spin. In the following, we discuss how the bath 
temperature and the strength of the applied magnetic 
field affect the decay of the elements of the reduced den- 
sity matrix. 
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FIG. 4: The decay of the componentsAa (t) and Ai(t) for 
g = (dotted lines), g = 5 (dashed lines), and g — 10 (solid 
lines). Other parameters are /3 = 1, 7 = A = 0, /i = 0.1a, 
A3(0) = i and Ai,2(0) = |. 

Clearly, if /i = 0, the vector component Xsit) vanishes 
when t —^ 00 regardless of the bath temperature. This 
means that /9ii(oo) = ^22(00) = ^, which is obviously in- 
dependent of the initial state of the central system. On 
the contrary, the off-diagonal elements tend asymptot- 
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ically to half of their initial values. This follows from 
the fact that the temperature-dependent quantity 7]°° is 
proportional to the magnetic field strength. The latter 
results are mainly due to the rotational symmetry of the 
model Hamiltonian together with the randomness of the 
interactions within the bath. Fig |4] illustrates the dif- 
ference in the decoherence process between the case of a 
static bath {g = 0) and a dynamic bath (g 7^ 0). We can 
see that the asymptotic value of Ai {t) decreases with the 
increase of g in contrast to X^lt) which assumes larger 
asymptotic values when g increases. However, at short 
times the above components decay slower with the in- 
crease of g, implying that strong quantum correlations 
within the environment suppress the effect of the deco- 
herence processes.. As we shall see below, the decay of 
the reduced density matrix exhibits a reverse behavior 
with respect to the temperature of the bath . This can 
be explained by the dependence of the decoherence time 
constant of our model, which turns out to be equal to 

r = Y^^if^ as revealed by Eqs. ^ and on 

the product gf3. Clearly, r ^ cxd as g ^ 00 or/and 
T 



T ^ 00 as g 
0, which confirms the above statements. 
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FIG. 5: Dependence of X3{t) and M{t) on the strength of the 
magnetic field in the case 00: fi = (dotted lines), /i — 

0.5a (dashed lines), and fi — 2a (solid lines). The parameters 
are 7 = 0, = 2, A3(0) = ^ and Ai,2(0) = |. 

Figure [5] illustrates the dependence of the components 
of the Bloch vector on the strength of the magnetic field. 
We can see that Xi{t) decays with the increase of fi 
whereas X^lt) approaches its initial value. Indeed, by 
making use of the following asymptotic expression of the 



incomplete gamma function2^ 

a- I (a-l)(a-2) 



r(a,z) - z-^-^e-^ 



1- 



z z" 
when z ^ 00 in I argz| < 37r/2, we obtain 

lim 7y°° = 1. 

fi,f3 — ^oo 



(91) 



(92) 



Therefore, if the ratio /i/a is infinitely big then the off- 
diagonal elements of the reduced density matrix tend 
asymptotically to zero; the diagonal ones assume their 
initial values. The above results can be explained by the 
fact that the effect of the bath on the dynamics of the 
central spin can be neglected when ^ is very large com- 
pared to a. The evolution in time is thus governed by the 
free Hamiltonian Hs which does not affect the diagonal 
elements of the reduced density matrix. The off-diagonal 
elements, however, show periodic oscillations; the vanish- 
ing asymptotic values obtained from Eqs. ([57)1 and ([^^ 
will never be reached since the decoherence time constant 
is infinite (a — > 0). 




-P 




0.4 
0.2 
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FIG. 6: Dependence of Xs(t) and \i{t) on the bath tempera- 
ture in the case N — > 00: f3 — (dotted lines), /3 = 1 (dashed 
lines), and /3 = 10 (solid lines). The parameters are 7 = 0, 
g = 2, fi = 0.5a, A3(0) = i and Ai,2(0) = |. 

From Fig. [6] it can be seen that the bath temperature 
has a reverse effect on the decay of the reduced density 
matrix elements. The components X3{t) and Xi{t) decay 
faster with the increase of T. Furthermore, we can see 
that the diagonal elements assume larger asymptotic val- 
ues in contrast with the off-diagonal ones. In the limit of 
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zero temperature, the asymptotic behavior is identical to 
the one corresponding to — > oo, see Eq. Indeed, 
at zero temperature the bath and the central spin evolve 
independently from each other as we already mentioned 
in the previous section. Once again, we find that the dy- 
namics of the central spin is governed by the free Hamil- 
tonian Hs which preserves the coherence of the central 
system. 




FIG. 7: The short-time behavior of \s,{t) and \\{t) in the case 
oo. The solid lines correspond to the exact solutions, 
the dotted lines denote the approximations (193194^ . Here, 
7 = 0, 3/3 = 1, /i = 0.5a, A3(0) = f and Ai,2(0) = |. 



Let us now discuss the short-time behavior of the re- 
duced dynamics. The aim here is to find simple analytical 
expressions which describe the variation of the reduced 
density matrix at short time scales. It is clear from the 
expressions of the functions ri{t), ({t) and ^{t) that the 
term of interest which describes the decay of the Bloch 

vector components is given by e (911+2) . We shall look 
for functions of the form e (3/5+2) g(^i^ where g{t) is some 
complex- valued function of the time . In the case of the 
off-diagonal elements, the ansatz g{t) — e^*^* can be jus- 
tified by the competition of two processes, namely oscil- 
lations due to the external magnetic field and damping 
due to the coupling with the environment. In the limiting 
case where the magnetic field is absent, it is found that 
for small values of the time, the decay is purely Gaussian. 
On the other hand if we assume that there is no coupling 
between the bath and the central spin, i.e. a = 0, then 
the dynamics is governed by the external magnetic field. 




FIG. 8: (Color online) Purity evolution for diflerent values 
of the bath temperature in the case N 00 with 7 = 0, 

g = 1 and fi — 0. The initial conditions are A3(0) — 

Ai,2(0) = i. 




FIG. 9: (Color online) Evolution in time of the purity for 
different values of the bath temperature in the case N 00 
with 7 = 0, g = 1 and jj. — a. 



The component A3 (t) is not affected by the magnetic field 
even when there is no coupling between the spin and the 
bath. It is shown in Ref. [la that this component decays 
two times faster than the other ones. Consequently, the 
short-time behavior of the reduced density matrix can be 
described by 



MO) 



exp 



V 2 + g(3J 



Ai(0)-zA2(0) 



exp 



.9/3 



(93) 
(94) 



In Fig. [21 the short time behavior of the Bloch vector 
components X^it) and Xi{t) is shown together with the 
approximations ([93|) and ([94|) : these are in good agree- 
ment with the exact solutions. 



13 




01234567 
at 



FIG. 10: Evolution of X^it) and \i{t) for iV = 600 (dotted 
lines) and N ^ oo (solid lines). The dashed lines correspond 
to the asymptotic values. Here, 7 = 2a, gr = 1, /3 = 1.5, 
A = 1, = 0.3q, A3(0) = i and Ai,2(0) = |. The plots cor- 
responding to Xait) are almost identical, the latter component 
saturates with respect to A*' faster than Ai(t). 

There exist many measures that allow for the quan- 
tification of the degree of the decoherence due to the in- 
teraction with an environment. In this work we use the 
measure D{t) = 1 — P{t), where 

P{t) = tr{p{tf} (95) 



is the purity of the central system. Note that in the pre- 
vious expression the trace is performed over the degrees 
of freedom of the central spin. The purity takes its max- 
imum value 1 at pure states; its minimum value, 1/2, 
corresponds to the fully mixed state p = I2/2. In our 
case, the purity can be expressed in terms of the Bloch 
vector components as 

P{t) = l{l + Ai + Ht? + AaW'). (96) 

The above expression shows that the decay of the pu- 
rity in the short-time regime described by Eqs. (|93p - ([94|) 
is Gaussian which reflects the non-Markovian character 
of the dynamics. The decay process is slowed down by 
decreasing the temperature of the bath and/or apply- 
ing a magnetic field of sufficient strength as illustrated 
in Figs. [5] and [51 When t 00, the central spin shows 
partial decoherence; if p — 0, then the asymptotic value 
of the purity is independent of the bath temperature as 
expected (see Fig. IH). 



B. The case 7 / O 

The time dependence of the Bloch vector components 
when the constant 7 is different from zero can be ob- 
tained with the same method used in the previous subsec- 
tion. Since 7 7^ 0, the quantities Mi^2 are m-dependent 
which means that the effect of the anisotropy constant 
has to be taken into account. When 7 7^ 0, we need to 
perform double integration with respect to the real vari- 
ables r and m as shown in Eq. ([7^ . By making use of 
the Riemann-Lebesgue lemma, it is possible to find the 
following asymptotic expression for the function ri{t) ob- 
tained by replacing p hy fi + 7m in Eq. ((78l) (see Fig. [TO)) 



00 

lim 77(t) = 1- ^{2 + gP)J2 + gPA f {p + -fmfe'^^'+''"'^'^3l3+2)-{2+gf!A)m'-^fQ^(^^^^^^2(^2 + gl3))dm. (97) 



Obviously, the functions ^(t) and (,{t) tend to zero 
when i — > 00. Hence, even if we set p = 0, the asymp- 
totic state is still temperature dependent. Nevertheless, 
the dependence of the Bloch vector components on the 
bath temperature is quiet similar to the one correspond- 
ing to 7 = 0. The influence of the magnetic field on the 
dynamics of the central spin is appreciable only when its 
strength is sufficiently large; this can be seen from the 
absence of oscillations in the components Xi(t) and Xsit) 



displayed in figure [Till Fig. [TT] shows that the off-diagonal 
elements decay slower and assume larger asymptotic val- 
ues when the anisotropy constant A increases. The op- 
posite situation holds for the component A3(i), that is 
when A decreases the latter component assumes larger 
asymptotic limits. 

At short times the diagonal elements of the reduced 
density matrix do not depend on A in contrast with the 
off-diagonal ones. In the case of the Heisenberg XY 
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FIG. 11: Dependence of X3{t) and \i{t) on the anisotropy constant in the case A'^ oo: A = (solid lines), A = 5 (dashed 
lines), and A = 10 (dotted lines). Here, y = 2a, gf3 ^ 1, fi ^ 0.1a, A3(0) = | and Ai,2(0) = |. 




FIG. 12: The short-time behavior of Xsit) and Xi{t) in the 
case N ^ oo. The solid lines correspond to the exact solu- 
tions, the the dotted lines denote the approximations (|98I99I) . 
The parameters are 7 = la, A = 0, gl3 — 1 and /i = 5a. The 
initial conditions are the same as in Fig llll 



model, i.e. when A = 0, the short-time behavior of 
the Bloch vector components can be determined with the 
same procedure used in the case where 7 = 0. The main 
difference here is that the contribution of the interac- 
tion V ~ -^=S^Jz has to be taken into account. Using 

the result we obtained in Eq. (|68p. we can describe the 
short-time behavior of the reduced density matrix by (see 



Fig.inD 



A3(0) 

l\2{t) 



Ai(o)-a2(o) 



cxp 



cxp 



\ 2 + gp) 



24-2 



2 + 5/3 



2iiit 



(98) 
(99) 



For A 7^ 0, the situation is much more complicated; here 
we only discuss the special case where /x = a = 0, and 
A >> 1. The last condition implies that the transverse 
term of Hb can be neglected compared to the longi- 
tudinal one. Under the above assumption, Hb simpli- 
fies to gA/NSl and thus all interactions are of Ising 
type. Therefore, the operators HgB and Hb commute 
with each other which means that the diagonal elements 
are not affected by the coupling to the environment. 
The coherence of the central spin can be calculated as 
usual. Taking the limit of an infinite number of spins 
and using the probability density function correspond- 
ing to the random variable to, we find that the off- 
diagonal elements decay according to the Gaussian law 

Hence the larger the anisotropy constant 



exp 1^ 2-1-3/3 A 

the slower the decay of the off-diagonal elements, which 
explains the behavior at short times of Ai {t) displayed in 
Fig. [nj More details about the case of Ising couplings 
can be found in Ref. H^- To end our discussion about 
the short-time behavior, it should be noted that the de- 
viation of the short-time expressions ([94|) and ((99|) from 
the exact solutions depends on the value of the strength 
of the magnetic field. For small values of /i, the above 
relations are valid at relatively large intervals of time. 
However, as /x increases, the domains of time for which 
the above approximations are valid become shorter. 

The variation in time of the purity in this case differs 
from the one corresponding to 7 = by the suppression 
of the damped oscillations caused by the external mag- 
netic field as shown in Fig [T31 This is mainly due to 
the interaction described by the Hamiltonian V . Conse- 
quently, the central spin decoheres less when 7 is equal to 
zero. The above result was expected because the longitu- 
dinal coupling vanishes: the central spin is less correlated 
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mainly due to the non-Markovian nature of the dynam- 
ics, which in turn follows from the time-independence of 
the bath correlation functions and the symmetry of the 
bath Hamiltonian. Also, it has be shown that the effect 
of low bath temperatures on the decoherence of the cen- 
tral spin is similar to that of strongly applied magnetic 
fields and large bath anisotropy . The results obtained 
in this work are valid for any number of spins in the 
environment and arbitrary values of the strength of the 
external magnetic field and the bath temperature. They 
are in good agreement with those of Ref. Hi] where the 
authors studied decoherence of electron spins in quantum 
dots. The model can be generalized to the case of two 
or more interacting qubits where questions related to the 
decoherence and the entanglement can be investigated. 



FIG. 13: (Color online) Purity evolution for different values of 
the bath temperature in the case N ^ oo with 7 = 2a, A = 1, 

g = 1, and /j, = a. The initial conditions are A3(0) — 

Ai.2(0) = i. 



to the environment and thus the destructive effect of the 
environment on the coherence of the two-level system 
is less appreciable. Indeed, the decoherence time con- 
stant is found to be inversely proportional to 7, namely 

„ 1 



4+2g;3 
2-1-7^(2+9/3)- 



This simply implies that r ^ as 
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CONCLUSION 



APPENDIX: DERIVATION OF THE 
ANALYTICAL FORM OF ri{t) 



In conclusion we have investigated the dynamics of a 
spin-i particle, subjected to the effect of a locally ap- 
plied external magnetic field, and coupled to anisotropic 
Heisenberg spin environment in thermal equilibrium. 
The reduced density matrix was analytically derived for 
finite number of spins in the environment and arbitrary 
values of the interaction strengths. The evolution in time 
of the central spin depends on the nature of interactions 
within the bath. In the case of ferromagnetic environ- 
ment, the decay of the Bloch vector components is Gaus- 
sian accompanied by fast damped oscillations. In the 
antiferromagnetic case, the components of the bloch vec- 
tor saturate with respect to the number of environmental 
spins and display partial decoherence. We showed that 
the partial trace over the degrees of freedom of the bath 
can be calculated using the convergence of the rescaled 
bath operators to normal independent Gaussian random 
variables. This allowed us to study the case of an infinite 
number of environmental spins, and to analytically derive 
the asymptotic behavior of the components of the Bloch 
vector. The above limit represents a good approxima- 
tion for the cases with finite number of spins {N ~ 100). 
At short time scales, the decay of the off-diagonal ele- 
ments is found to be Gaussian with a decoherence time 

This result is 



constant given by r 



2+g/? 



(7 = 0). 



This appendix is devoted to the derivation of the 
asymptotic behavior and the analytical form of the func- 
tion ri{t) appearing in Eq. ((79| . Explicitly we have 
(A-0) 



vit) 



Z 



g-(3/3+2)r^ 



■ sm' 



(A.l) 

By making the following change of variable — — 
fi^ and taking into account the trigonometric equality 



. 2 

sm X - 
tion as 



vit) 



1 — cos(2a;) 



we can rewrite the above func- 



dv 



1 — cos(2wi) 
t^^-i2+g0)v- (1 - cos(2i;i))] }• (A.2) 



The Riemann-Lebesgue lemma implies that the second 
and the fourth terms involving the cosine function in the 
above expression vanish when t ^ 00. The first term can 
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be easily evaluated and we simply get 



2 + g(3 



where we have made the change of variable (2 + gP)v'^ 

. Taking into account the expression of Z in Eq. ([7^ 
we obtain the asymptotic expression of •q{t) displayed in 
(A.3) Eq. 



The third term reads 

oo 

dv 



2 / Z-Q-i'^+aP)^' ^ 
J V J V 



^e-"^=r(0,(2+g/3V^ 



(A. 4) The second term simplifies to 



Re |2 e j v dv exp 



-(2 + 3/3) « 



it 



2 + gP 



)']} 



e ^+sPRe< / dv 



2it 



2 + gP {2 + g(3)i 

6^ S 



dv e 



(A.5) 



l2 



where 6 = \/2 + .g/3(/i - 



2+9/3' 



One can easily check that denominator of the quantity under the sign of integral 



Re (/i) 



1 



2 + 5/3 



exp 



-{2^gl3)p? 



.9/3 



cos^2/it 



(A.6) 

The second integral is given by the complementary error 
function, namely 



{2 + g^i 



-erfc 



(/.+ ^)y2T^]. (A.7) 



It is then sufficient to use the property 2 Im erfc(a + ii) = 

i erf(a+ii) — erf (a— it) , where a is real and Im(a;) stands 

for the imaginary part of x, to get the first three terms 
appearing in the right-hand side of Eq. ([7^ . 
Similarly, we have 



Re<^ 2 



dv 



■ exp 



-{2 + g[3){v 



Re{ 2 



it 



2 

ds 



5/3) ]} 

}, (A. 



where we have introduced the new variable s — (w + 
)-v/2 + g(3- By multiplying the numerator and the 



\/2q^ 



given by 



we get two new integrals. The first one is 



Re<^ 2 



ds- 



e 2+3(3 Re 



2+9/3 



^e-| 
s i 



(A.9) 



where ^2 = (2 + gP)ti^ + 2^it. The remaining integral 
defines the function Al, namely 



M{t;fi,f3) = exp| 



OQ o 

— — dsV (A.IO) 



The analytical expressions of the functions ^(t) and ({t) 
can be determined with the same method. In the case 
7 ^ we should replace /i by /i + 7m and then perform 
the integration with respect to m. For practical investi- 
gation, numerical integration is used. 
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